arXiv:1502.07544vl [cond-mat.supr-con] 26 Feb 2015 


Strong Interplay between Stripe Spin Fluctuations, Nematicity and Superconductivity in FeSe 


Qisi Wang , 1 Yao Shen , 1 Bingying Pan , 1 Yiqing Hao , 1 Mingwei Ma , 2 Fang Zhou , 2 
P. Steffens , 3 K. Schmalzl , 4 T. R. Forrest , 5 M. Abdel-Hafiez , 6,7 D. A. Chareev , 8 
A. N. Vasiliev, 9 - ia 11 P Bourges , 12 Y. Sidis , 12 Huibo Cao , 13 and Jun Zhao * 1,14 

1 State Key Laboratory of Surface Physics and Department of Physics, Fudan University, Shanghai 200433, China 
2 Beijing National Laboratory for Condensed Matter Physics, 

Institute of Physics, Chinese Academy of Science, Beijing 100190, China 
3 Institut Laue-Langevin, 71 Avenue des Martyrs, 38042 Grenoble Cedex 9, France 
4 Juelich Centre for Neutron Science JCNS Forschungszentrum Juelich GmbH, Outstation at ILL, 38042 Grenoble, France 
5 European Synchrotron Radiation Facility, BP 220, F-38043 Grenoble Cedex, France 
6 Center for High Pressure Science and Technology Advanced Research, Shanghai, 201203, China 
7 Faculty of Science, Physics Department, Fayoum University, 63514 Fayoum, Egypt 
8 Institute of Experimental Mineralogy, Russian Academy of Sciences, 142432 Chernogolovka, Moscow District, Russia 
9 Low Temperature Physics and Superconductivity Department, 

M.V. Lomonosov Moscow State University, 119991 Moscow, Russia 
10 Theoretical Physics and Applied Mathematics Department, 

Ural Federal University, 620002 Ekaterinburg, Russia 

11 National University of Science and Technology “MISiS ”, Moscow 119049, Russia 

12 Laboratoire Leon Brillouin, CEA-CNRS, CEA-Saclay, 91191 Gif sur Yvette, France 

13 Neutron Scattering Science Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6393, USA 
14 Collaborative Innovation Center of Advanced Microstructures, Fudan University, Shanghai 200433, China 


Elucidating the microscopic origin of nematic order in iron-based superconducting materials is 
important because the interactions that drive nematic order may also mediate the Cooper pairing 1 ! 
Nematic order breaks fourfold rotational symmetry in the iron plane, which is believed to be driven 
by either orbital or spin degrees of freedom 15 . However, as the nematic phase often develops at 
a temperature just above or coincides with a stripe magnetic phase transition, experimentally de¬ 
termining the dominant driving force of nematic order is difficult 1 Here, we use neutron scat- 
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tering to study structurally the simplest iron-based superconductor FeSe (ref. [7]), which displays a 
nematic (orthorhombic) phase transition at T s = 90 K, but does not order antiferromagnetically. 
Our data reveal substantial stripe spin fluctuations, which are coupled with orthorhombicity and 
are enhanced abruptly on cooling to below T s . Moreover, a sharp spin resonance develops in the 
superconducting state, whose energy (~ 4 meV) is consistent with an electron boson coupling mode 
revealed by scanning tunneling spectroscopy®, thereby suggesting a spin fluctuation-mediated sign¬ 
changing pairing symmetry. By normalizing the dynamic susceptibility into absolute units, we show 
that the magnetic spectral weight in FeSe is comparable to that of the iron arsenides 910 . Our findings 
support recent theoretical proposals that both nematicity and superconductivity are driven by spin 
fluctuation^^®!. 

Most parent compounds of iron-based superconductors exhibit a stripe-type long-range antiferromagnetic 
(AFM) order which is pre-empted by a nematic order: a correlation of electronic states which breaks 
rotational, but not translational, symmetry. Superconductivity emerges when the magnetic and nematic 
order are partially or completely suppressed by chemical doping or by the application of pressure®. The 
stripe AFM order consists of columns of parallel spins along the orthorhombic b direction, together with 
antiparallel spins along the a direction. Similar to the stripe AFM order, the nematic order also breaks the 
fourfold rotational symmetry, which is signaled by the tetragonal to orthorhombic structure phase transition 
and pronounced in-plane anisotropy of electronic and magnetic properties® 15 18 . It has been proposed that 
nematicity could be driven either by orbital or spin fluctuations, and that orbital fluctuations tend to lead 
to a sign-preserving ,v ++ -wave pairing, while spin fluctuations favor a sign-changing s^-wave or d -wave 
pairing 1 6 14 19 20 . However, as orbital and spin degrees of freedom are coupled and could be easily affected 
by the nearby stripe magnetic order, it remains elusive which of them is the primary driving force of 
nematicity 151419 . 

FeSe (T c « 8 K) has attracted great attention not only because of the simple crystal structure (Fig. la), 










3 


but also because it displays a variety of exotic properties unprecedented for other iron based superconduc¬ 
tors. For example, the T c of FeSe increases to ~40 K under pressure 21 or by ion/molecule intercalation 22 . 
In addition, the T c of single layer FeSe thin film is as high as 100 K, which is significantly higher than in 
other iron based superconductors 23 . More interestingly, unlike most iron-based materials, the tetragonal 
to orthorhombic structural transition in bulk FeSe is not followed by a stripe magnetic order 7 , providing 
an exciting opportunity to elucidate the microscopic origin of nematicity and its interplay with supercon¬ 
ductivity. The absence of stripe magnetic order in FeSe seems to cast doubt on the spin driven nematicity 
scenario. Moreover, recent nuclear magnetic resonance (NMR) measurements suggested that there were 
little spin fluctuations above T s in the tetragonal phase, which was also interpreted as a breakdown of the 
spin scenario 24 25 . However, NMR only probes momentum-integrated spin fluctuations at very low energies 
(~ O.lyueV or lower). The momentum dependence of the higher energy spin fluctuations-especially at the 
energy scale close to the superconducting gap, which is believed to be more important in driving nematicity 
and superconductivity 1114 -remains unknown. This issue could be addressed by inelastic neutron scattering 
measurements that probe spin fluctuations over a wide range of momentum and energy. 

Neutron scattering studies on FeSe single crystals have been hampered by the lack of high quality sam¬ 
ples with the correct phase. Recently, advances in crystal growth with vapor-transport and floating zone 
techniques have allowed us to grow FeSe single crystals which are significantly larger than what was pre¬ 
viously available 26 27 . The superconducting properties of our sample were characterized by DC magnetic 
susceptibility and resistivity measurements which show an onset T c of 8.7 K with a transition width of ~0.3 
K, indicating the high quality of the sample (Figs, lb, lc). Clear kinks on magnetic susceptibility and 
resistivity associated with the tetragonal to orthorhombic structure transition are also observed close to 90 
K. 

We first use elastic neutron scattering to study the structural and magnetic ordering properties of our FeSe 
samples. A broadening of the (4, 0, 0) structural peak is observed below 90 K, indicative of the structural 
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phase transition from the tetragonal to orthorhombic symmetry (Fig. 2a). The broadened line shape can 
be fitted with two Gaussian peaks since the sample has orthogonal twin domains and both (4, 0, 0) and 
(0, 4, 0) peaks are covered by the scan. No significant change of the peak width is observed across T c 
within our instrumental resolution. On the other hand, no magnetic Bragg peaks associated with the stripe 
or double stripe magnetic order are observed (not shown) at temperatures down to 1.5 K, consistent with 
previous measurements of powder samples 7 . Instead, in the inelastic channel, we have observed strong 
spin fluctuations near (1, 0, 0), which corresponds to the stripe AFM wavevector of the parent compounds 
of iron-based superconductors 6 ! To determine the momentum dependence of the spin fluctuations and 
their interplay with superconductivity, we performed rocking/transverse and radial/longitudinal (the scan 
directions are perpendicular and along Q, respectively) Q-scans below and above T c . As shown in Figs. 
2b, 2d, representative Q-scans at 4 meV are commensurate near (1, 0, 0) at T = 11 K in both transverse and 
longitudinal directions with no observable anisotropy. The peak intensity is drastically enhanced below T c , 
which is reminiscent of a magnetic resonant mode observed in other iron-based superconductor^ 9 10 14 28 29 i 
Conversely, the scattering at 2.5 meV is suppressed upon entering the superconducting state due to the 
opening of the superconducting spin gap (Fig. 2c). The redistribution of the magnetic spectral weight 
across T c indicates that the spin excitations near (1, 0, 0) are closely related to superconductivity. In order 
to clarify the detailed momentum structure of the superconductivity-induced magnetic excitations, we have 
subtracted the signal of the normal state from that of the superconducting state and plot a 2D contour map 
interpolated from a series of Q-scans at 4 meV (Fig. 2e). The outcome shows that the spin excitation 
spectra are very sharp with little anisotropy (within our instrumental accuracy). In addition to the results 
shown near (1,0, 0), we also performed similar measurements in the second magnetic Brillouin zone near 
(2, 1, 0) associated with the stripe magnetic structure (Fig. 2e). A similar signal is also observed, but with 
weaker intensity because of the decreased magnetic form factor. These results unambiguously demonstrate 
that the scattering that we observe here is pure magnetic fluctuation associated with stripe magnetism rather 
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than phonons as the scattering strength from phonons is related to (Q-£) 2 , where £, is the polarization vector 
of the phonon. 

Figure 3 summarized the energy dependence of dynamic spin correlation function S (Q ,co) at Q=(l, 0, 0) 
at different temperatures. The figure confirms that the spectral weight loss in the superconducting spin gap 
(< 3 meV) is compensated by a sharp resonance mode at around 4 meV. Moreover, the detailed temperature 
dependence of the scattering at 4 meV shows an order-parameter-like behavior and is clearly coupled to 
the occurrence of superconductivity (Fig. 4a). The spin resonance mode has been interpreted either as a 
spin exciton within the superconducting gap arising from scattering between portions of the Fermi surface 
where the superconducting gap function has an opposite sigrf 14 ^ or as a broad hump structure induced by 
overshoot in the magnetic spectrum above the superconducting gap in a sign-preserving s ++ pairing statcf 2f) [ 
The sharp mode that we observed here is consistent with the spin exciton model as the mode energy (4 meV) 
is below the superconducting gap (2A«5 meV) (ref. 1301). and the energy width (~ 1.2 meV) of the mode is 
essentially resolution-limited and much sharper than in other iron-based superconductor^ 9 10 142829 . More 
interestingly, the resonance energy (E r =4meY^5.3k B T c ) is consistent with the electron boson coupling 
mode (~3.8 meV) revealed by scanning tunneling spectroscopy 8 , thereby suggesting strong electron-spin 
excitations coupling in this system. These results are consistent with a spin fluctuation-mediated sign 
changing pairing mechanism, but inconsistent with an orbital fluctuations-mediated sign-preserving ,v ++ - 
wave pairing mechanism! 4 20 . 

Although sharp and commensurate stripe spin fluctuations persist at all temperatures measured, the sys¬ 
tem remains paramagnetic at low temperature. Theoretically, it has been shown that the magnetic interac¬ 
tions in FeSe are much more frustrated than in iron arsenides and therefore prevent long-range magnetic 
order 1112 . Hence, it is informative to compare the magnetic spectral weight in FeSe with that of iron 
arsenide superconductors. We have calculated absolute units of the imaginary part of the dynamic sus¬ 
ceptibility x"(Q,u>) by normalizing S(Q.co) for the thermal population factor and the intensity of acoustic 
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phonons (Fig. 3b and Supplementary Information). The outcome reveals that the integrated resonance spec¬ 
tral weight (0.00212 /r|/Fe) is about 30% of that in the carrier doped BaFe 1 . 8 5 Coo.i 5 As 2 (E r = 9.5 meV) 
(ref. 0), but two times larger than the damped resonance mode in the isovalently doped BaFc| 85 Ru 0 i S As 2 
(T c =14 K, E r - 5.5 meV) (rcf.lTOl). Since the T c (8.7 K) of FeSe is also about a factor of three lower than in 
BaFei. 85 Co 0 15 AS 2 (T c = 25 K), the overall magnetic spectral weight in both systems should be comparable. 

Having established the interplay between the spin fluctuations and superconductivity, we now turn to 
the impact of nematicity on the spin fluctuations. Previous NMR measurements suggested the absence of 
spin fluctuations above T s in the tetragonal phase 2425 . In contrast, our neutron scattering measurements 
show substantial spin fluctuations in the tetragonal phase (T = 110 K) (Figs. 3a, 2d). We note that the 
energy dependence of the dynamical spin correlation function S (Q ,of) displays a spin gap-like feature at 
low energies at T = 110 K (Fig. 3a), which is confirmed by the featureless Q-scan at 2.5 meV (Fig. 2c). 
These results agree with a theoretically predicted gapped nematic quantum paramagnetic state with low 
carrier density in FeSe (ref. ITTl) . which naturally accounts for the absence of low energy spin fluctuations 
above T s suggested by NMR measurements 2425 . The most striking observation is that the spin fluctuations 
are enhanced abruptly in the orthorhombic phase at T = 11 K (Fig. 3a). We note that the increase 
of the spin fluctuations is more pronounced at low energies. To determine if the increase of the spin 
fluctuation is indeed associated with the nematic order, we carefully measured the temperature dependence 
of the scattering at 2.5 meV, which is the lowest energy that can be measured in our thermal triple axis 
spectrometer with a reasonable background. Intriguingly, a comparison of the temperature evolution of 
the S(Q,a>) with the orthorhombicity S(T) = (a - b)/(a + b ) reveals that the enhancement of the S(Q,a») 
is clearly coupled to the development of the nematic (orthorhombic) phase (Fig. 4b). These results are 
consistent with the recent proposals (based on either itinerant or local moment pictures) that the nematic 
order is driven by spin fluctuations 1 2 11 13 . In a local moment model that frustrated magnetic interactions 
drive nematic order in FeSe, once the orthorhombic distortion develops, the effective nearest-neighbor 
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exchange couplings J lx and J ly become non-equal and the frustration is partially released, therefore making 
the system move toward the stripe ordered phase. As a result the spin fluctuations at the stripe ordering 
wavevector are enhanced. 

It is interesting to compare the spin fluctuations of FeSe with that of iron selenide superconduc¬ 
tors without nematic order. The low energy spin fluctuations in FcTci_ a Sc a and Rb A Fc 2 - v Sc 2 appear at 
Q=(l, -0.3 < £ < 0.3) and Q=(l, +0.5), respectively 28 29 31 . Different from FeSe, the dynamic spin corre¬ 
lation S (Q.oj) of FeTC|_ A Se A displays little temperature dependence from T c to 300 K (ref. [28]). Moreover, 
the spin fluctuations of FeTei_ Y Se A are broad and incommcnsuratc/anisotropid- 28 29 , in contrast to the sharp 
and commensurate spin fluctuations at the stripe AFM wavevector in FeSe. Therefore, FeSe is closer to the 
stripe magnetic instability and consequently with a larger spin-spin correlation length. These results further 
imply that nematicity is driven by stripe spin fluctuations, though superconductivity can be mediated by 
spin fluctuations either at or away from the stripe AFM wavevector. 

In summary, we have reported evidence of strong coupling between the stripe spin fluctuations, ne¬ 
maticity and superconductivity in single crystalline FeSe. Contrary to earlier NMR measurements 2425 , our 
neutron scattering data reveal substantial commensurate stripe spin fluctuations in the tetragonal phase, 
which are further enhanced in the nematic phase. Moreover, a resolution-limited sharp spin resonance ap¬ 
pears well below the superconducting gap and is coupled with electronic density of states, indicating a spin 
fluctuations-mediated sign-changing pairing symmetry rather than an orbital fluctuations-mediated sign¬ 
preserving s ++ -wave pairing symmetry. These results are in agreement with the theoretical predictions that 
nematicity and superconductivity are driven by spin fluctuations 1 2 11 14 . We believe that the elucidation of 
the interplay between spin fluctuations, nematicity and superconductivity will have important implications 
for the understanding of other exotic properties of FeSe, such as the drastically increased T c under external 
pressure or substrate strain 21 23 . 

Note added: After we finished this paper, we became aware of a related preprint describing neutron 
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scattering measurements on FeSe powder samples 32 . 

* Correspondence and requests for materials should be addressed to J.Z. (zhaoj@fudan.edu.cn). 
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Figure 1: Orthorhombic crystal structure, magnetic susceptibility and resistivity of FeSe single crystal, a, Schematic diagram of FeSe 
crystal structure, b. The DC magnetic susceptibility measurements on the single-crystalline FeSe sample. A sharp superconducting transition 
is observed at T c = 8.7 K in the ZFC measurement in a magnetic field of H = 10 Oe, indicating ~ 100% exclusion of the external magnetic 
field. The screening is slightly larger than -1 because of the demagnetization effect. The inset shows the susceptibility measured in a magnetic 
field of H = 20 kOe. The magnetic fields are applied perpendicular to the c axis, c, In-plane resistivity as a function of temperature. The inset 
shows data around T s = 90 K on an enlarged scale. 
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Figure 2: Structure phase transition and momentum dependence of the spin fluctuations at various temperatures in FeSe. The inelastic 
neutron scattering measurements were carried out on the IN20 thermal triple axis spectrometer at the Institut Laue-Langevin, Grenoble, France 
and the 2T1 thermal triple axis spectrometer at the Laboratoire Leon Brillouin, France. The FeSe single crystals are co-aligned in the (H, K, 0) 
horizontal scattering plane within ~ 3 degrees mosaicity for the measurements. The elastic measurements were performed on one piece of small 
single crystal on the HB-3A four-circle single-crystal diffractometer at the High-Flux Isotope Reactor at the Oak Ridge National Laboratory, 
United States (The instrument configurations are described in the Supplementary Information). We present the data by defining the wave vector 
Q at (q x , q y , q.) as (h, k, /)= (q x ci/2n, q y a/2n, q z c/2n) reciprocal lattice units (r.l.u.) in the orthorhombic unit cell, a, Temperature dependence 
of the (4, 0, 0)/(0, 4, 0) nuclear reflections. The Bragg peak is significantly broadened below the tetragonal-to-orthorhombic phase transition, 
b-d, Q-scans near (1, 0, 0) at various energies and temperatures; linear backgrounds are subtracted (see the Supplementary Information). The 
scan directions are marked by green arrows in the insets, e, 2D contour plot of the temperature difference scattering [5(1.5^) - S(llA')] 
interpolated from a series of Q-scans at 4 meV. The error bars indicate one standard deviation. 
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Figure 3: Energy dependence of spin excitations for FeSe in the superconducting state (T = 1.5 K) and normal state ('/' = 11 and 110 

K) a, Energy dependence of the dynamic spin correlation function S (Q,cu) at Q =(1, 0, 0) after a background correction. The background is 
measured at Q = (0.944, 0.330, 0) and Q = (0.944, -0.330, 0), 0) (see the Supplementary Information). The open circles are data fitted with 
Q-scans. b, Energy dependence of the imaginary part of the dynamic susceptibility ;f"(Q,tu) in the superconducting state (T = 1.5 K), and the 
normal state (T = 11 K). The data are obtained from 5(Q ,a>) by correcting for the Bose-population factor and are normalized to absolute units 
with acoustic phonons as described in the Supplementary Information. The solid curves are guides to the eye. The shaded area denotes the 
resonance spectral weight. The error bars indicate one standard deviation. 













14 




Figure 4: Temperature dependence of spin fluctuations in FeSe. a, Temperature dependence of dynamic spin correlation S (Q,ai) at E = 4 
meV, which clearly shows a kink at T c . The inset displays the temperature evolution of^'"(Q,w). b. Temperature dependence of S(Q,at) at 
E = 2.5 meV and the orthorhombicity <5(7’) = (a - b)/(a + b) shows an order-parameter-like behavior with an onset at T s . The orthorhombicity 
is adapted from the X-ray diffraction data in ref.[7] The inset shows the temperature evolution of^"(Q,w) which also displays a kink at T s . We 
note that the decrease of the scattering intensity at 1.5 K is due to the opening of the superconducting spin gap (Fig. 3). The open circles are 
data fitted with Q-scans. The solid lines are a guide to the eye. The error bars indicate one standard deviation. 
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Supplementary Information: Strong Interplay between Stripe Spin Fluctuations, Nematicity and 
Superconductivity in FeSe 

I. Instrument configurations for the elastic and inelastic neutron scattering experiments 

Elastic neutron diffraction on FeSe single crystals was measured on the HB-3A four-circle diffractometer 
at the High Flux Isotope Reactor at Oak Ridge National Faboratory, United States. The neutron energy 
of 34.4 meV was used from a bent perfect Si-220 monochromator 1 . We used the high resolution mode 
by flatting the monochromator (mbend=50) in order to track the lattice distortions (Fig. 2a). The inelas¬ 
tic neutron scattering measurements were carried out on the IN20 thermal triple axis spectrometer at the 
Institut Faue-Fangevin, Grenoble, France, and the 2T1 thermal triple axis spectrometer at the Faboratoire 
Feon Brillouin, France. For the measurements performed on IN20 (Figs. 2b-2e, 3a, 3b, and 4b), we used 
a focusing Si(lll) as monochromator and a pyrolytic graphite [PG(002)] as analyzer. This setup yielded 
an energy resolution of about 1 meV at (1, 0, 0) at E = 0 meV. For the measurements performed on 2T1 
(Fig. 3a), PG(002) was used as the monochromator and analyzer. A PG filter was installed in front of 
the analyzer to eliminate the contamination from the higher-order neutrons. A correction was also made 
for monitor contamination by higher-order neutrons. For both triple axis spectrometers, the final neutron 
energy was fixed at F/= 14.7 meV and no collimation was used. 

II. Raw data, background subtraction and absolute units normalization 

For typical inelastic neutron scattering experiments, the background is momentum, energy and tempera¬ 
ture dependent. SFig. 1 shows several representative raw Q-scans measured at various temperatures. Each 
Q-scan can be fitted by a single Gaussian peak on a linear background. The data presented in Fig. 2 are 
obtained by subtracting the linear background from the raw Q-scans. 

SFig. 2 shows several temperature difference Q-scans [5(1 ,5K) - 5(1 \K)\ at 4 meV near the first Bril¬ 
louin zone center (1, 0, 0) and the second Brillouin zone center (2, 1, 0). The 2D contour plot in Fig. 2e 
was interpolated from a series of such Q-scans. 
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SFig. 3 summarizes the raw energy scans at 1.5 K, 11 K and 110 K. The background was estimated as 
the average intensity at Q = (0.944, 0.330, 0) and Q = (0.944, -0.330, 0). The background-subtracted data 
are presented in Fig. 3. 

SFig. 4a shows the temperature dependence of the scattering at the signal [Q=(l, 0, 0)] and background 
positions at 2.5 meV. As expected, the scattering intensity of the background decreases gradually with 
decreasing temperature. Nevertheless, on cooling to below T s =90 K, a sudden increase of the scattering 
intensity at the signal is clearly seen (SFig. 4a). Similar behavior is also observed near T c = 8.7 K at 4 
meV (SFig. 4b). 

The absolute units of the imaginary part of the dynamic susceptibility _y"(Q,6u) (Fig. 3b) were calculated 
by comparing the intensity of spin fluctuations with that of the acoustic phonons. This approach has been 
used to normalize the spin fluctuation intensity in several iron based superconductors 2 3 and is intensively 
discussed in ref. 4 

* Correspondence and requests for materials should be addressed to J.Z. (zhaoj@fudan.edu.cn). 
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SFig 1: Representative raw Q-scans measured at various temperatures in FeSe. The Q-scan can be fitted by a single Gaussian peak on 
a linear background. The dashed line indicates the background, a, Rocking scan at E=4 meV, Q=(l, 0, 0). T=1.5 K. b. Rocking scan at E=4 
meV, Q=(l, 0, 0), T=ll K. c, Rocking scan at E=4 meV, Q=(l, 0, 0), T=110 K. d. Rocking scan at E=2.5 meV, Q=(l, 0, 0), T=ll K. e, 
Rocking scan at E=2.5 meV, Q=(l, 0, 0), T=110 K. 
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SFig 2: Representative temperature difference Q-scans [S(1.5K) - S(11K)] at 4 meV near (1, 0, 0) and (2, 1, 0). The 2D contour plot 
in Fig. 2e was interpolated from a series of such Q-scans. The scan directions are marked in the insets. Each scan can be fitted by a single 
Gaussian peak. No significant anisotropy of the peak width is observed, a, Rocking scan at E=4 meV, Q=(l, 0, 0). b, Hscan at E=4 meV, 
Q=(l, 0, 0). c, Hscan at E=4 meV, Q=(2, 1, 0). d, Kscan at E=4 meV, Q=(2, 1, 0). 
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SFig 3: Energy dependence of the scattering at the signal [Q=(l, 0, 0)] and background positions. The background was estimated as 
the average intensity at Q = (0.944, 0.330, 0) and Q = (0.944, -0.330, 0). a, T=1.5 K. b, T=ll K. c, T=110 K. The overall magnetic spectral 
weight is clearly enhanced on cooling from llOKto 11 Kat the energies measured. 
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SFig 4: Temperature dependence of the scattering at the signal [Q=(l, 0, 0)] and background positions at 2.5 meV and 4 meV. The 

background-subtracted data are presented in Fig. 4. a, Temperature dependence of the scattering at the signal [Q=(l, 0, 0)] and background at 
2.5 meV. The background was estimated as the average intensity at Q = (0.944, 0.330, 0) and Q = (0.944, -0.330, 0). Although the background 
decreases gradually with decreasing temperature, the signal exhibits a sudden increase at T s =90 K. We note that the decrease of the scattering 
intensity at the signal at 1.5 K is simply due to the opening of the superconducting spin gap. b, Temperature dependent data for 4 meV with the 
background measured at Q=(l, 1, 0) and Q=(l, -0.6, 0). Since this scan was measured in a relatively narrow temperature range (3K to 21K), 
the background was estimated by a linear fitting of the data points collected at eight temperatures (blue squares). This is justified as most data 
points fall on the fitting curve (blue solid line). The data presented in SFig. 4b were collected on 2T1. All other inelastic neutron scattering 


data were collected on IN20. 





